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Abstract 

Recently an upper bound on the rate of the lepton number violating 
decay K + — > /x + /i + 7r~ has been significantly reduced by the E865 experi- 
ment at BNL and further improvement is expected in the near future. We 
study this process as a possible source of information on neutrino masses 
and mixings. We find that it is insensitive to the light (eV domain) and 
heavy(GeV domain) neutrinos. However due to the effect of a resonant 
enhancement this decay is very sensitive to neutrinos Vj in the mass region 
245 MeV < m u . < 389 MeV. At present experimental sensitivity we de- 
duce new stringent limits on the neutrino mixing matrix element U^j for 
neutrino masses in this region. 



1 Introduction 



Nowadays there are sufficient indications to believe that neutrinos are massive 
particles mixing with each other jl]. These indications come from both experi- 
mental and theoretical sides. The solar neutrinos deficit, the atmospheric neutrino 
anomaly and the results of the LSND neutrino oscillation experiment, all can be 
explained in terms of neutrino oscillations implying non-zero neutrino masses and 
mixings. On the theoretical side almost all phenomenologically viable models of 
the physics beyond the standard model(SM) predict non-zero masses for neutrinos 
which can be either Majorana or Dirac particles. 

1 On leave from the Joint Institute for Nuclear Research, Dubna, Russia 



1 



Majorana masses violate the conservation of total lepton number by two units 
AL = 2. Therefore lepton number violating^ ) processes represent a most ap- 
propriate tool to address the question of whether neutrinos are Majorana or Dirac 
particles. Various Ijj processes have been studied in the literature in this respect. 
Among them there are the neutrinoless nuclear double beta (0^/5/3) decay 0, §, 
the decay K + — > ||, || @, [7], §|, nuclear muon to positron |J or to an- 



tiunion JT0[ conversion, trimuon production in neutrino-nucleon scattering [I I 



the process e p — > ulil^X, relevant for HERA |12|], as well as direct produc- 
tion of heavy Majorana neutrinos at various colliders [13]. The analysis made 
in the literature || [b| leads to the conclusion that if these processes are me- 
diated by the Majorana neutrino exchange then, except 0z//?/3-decay, they can 
hardly be observed experimentally. This analysis relies on the current neutrino 
oscillation data, and on certain assumptions about the neutrino mass matrix. In 
the present paper we concentrate on the neutrinoless double muon decay of kaon 
K + — > fi + {i + 7i~. We will show that despite the above conclusion being true for 
contributions of the neutrino states much lighter or much heavier than the typical 
energy of the K + — > p, + fi + 7i~ decay, there is still a special window in the neutrino 
sector which can be efficiently probed by searching for this process. This window 
is in the neutrino mass range 245 MeV < m Vj < 389 MeV, where the s-channel 
neutrino contribution to the K + — > fi + fi + 7i~ decay is resonantly enhanced, there- 
fore making this decay very sensitive to the neutrinos in this mass domain. If 
neutrinos with masses in this region exist, then from present experimental data 
we can extract stringent limits on their mixing with v^. We derive these limits 
from the upper bound on the branching ratio of K + — > fi + fi + 7T~ recently obtained 
by E865 experiment at BNL [pl| . 

2 K + —> /i + /i + 7r~ decay in Standard Model with 
Majorana neutrinos 

In the SM extension with Majorana neutrinos there are two lowest order diagrams, 
shown in Fig.l, which contribute to K + — > fi + fi + 7i~ decay. These diagrams were 
first considered long ago in Refs. [§, |j. Here we are studying previously over- 
looked aspects of this decay. We concentrate on the s-channel neutrino exchange 
diagram in Fig. 1(a) which plays a central role in our analysis. The t-channel 
diagram in Fig. 1(b) requires in general a detailed hadronic structure calculation. 
In Ref. || this diagram was evaluated in the Bethe-Salpeter approach and shown 
to be an order of magnitude smaller than the diagram in Fig. 1(a), for light and 
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intermediate mass neutrinos. As we will see, in the neutrino mass domain of our 
main interest, the diagram in Fig. 1(a) absolutely dominates over the t-channel 
diagram in Fig. 1(b), independently of hadronic structure. 

The contribution from the factorizable s-channel diagram in Fig. 1(a) can be 
calculated in a straightforward way, without referring to any hadronic structure 
model. A final result for the K + — > /! + /i + 7r~ decay rate is given by 
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The unitary mixing matrix Uij relates v[ = UijVj weak v' and mass v neutrino 
eigenstates. The numerical constant in Eq. ([]]) is c = (G F / 32) (tt)~ 3 f 2 f 2 m 5 K \V u d\ 2 \ V u , 
where f K = 1.28 f^, f v = 0.668 m n and mx = 494 MeV is the K-meson mass. 
The functions G(z) and H(z±, z 2 ) in Eq. ([!]) after the phase space integration can 
be written in an explicit algebraic form 
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H(z 1 ,z 2 ) = h__(z 1 )h__(z 2 ) +xl[r + (z 1 z 2 ) - x 2 ^, z 2 , 1)] - r ^{zxz 2 )t{z u z 2 , 



Here we defined X; L = mi/m K and h±±(z) = z ± x 2 ± x 2 , r±{z\z 2 ) = z\z 2 — 
4 ± 4> t( Zl ,z 2 ,z 3 ) = Z! + z 2 - 2z 2 , <f>{z) = \ 1/2 (l,xl,z)\ 1 / 2 (z 1 x 2 l ,xl) with 
X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2xz. The integration limits in Eq. ([!]) 
are 
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with y = si/m 2 . 

First we assume that neutrinos can be separated into light v\. and heavy Nk 
states, with masses m vi « J~sf[ and \J~s{ « M^k- Then the Eq. ([[]) can be 
approximately written as 
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where the average neutrino masses are determined in the standard way 



KW= E tf*™*, (m^)^= E Kk M Nl- (5) 

k=light k=heavy 

The approximate formula (|J) can be used for extracting limits on the average 
neutrino masses from the experimental data. This leads to the limits: | (m u ) \ < 
Exp[y\ ((M^ 1 )! < Exp(N). We point out that in this case the upper bounds 

must satisfy the consistency conditions Exp(v) « \J~s~[ ~ m K , Exp(N)^ 1 » 

\J~sf ~ m K . Otherwise the so derived limits are not applicable to (m u ), (Mjy 1 ) 
as it is in Refs. [|7|, || If consistency conditions are not satisfied one has to 
use the initial formula ([!]). Similar consistency conditions take place for the other 
Jjj processes. 

The dimensionless coefficients in Eq. (^) are A v = 4.0 x 10~ 31 , An = 7.0 x 
10~ 32 , A u n = 1-7 x 10~ 31 . With these numbers we can estimate the current 
upper bound on the K + — > fi + fi + 7r~ decay rate from the experimental data on 
other processes. 

Atmospheric and solar neutrino oscillation data, combined with the tritium 
beta decay endpoint, allow one to set upper bounds on the masses of the known 
three neutrinos [1(J m e ^ tT < 3 eV. Thus in the three neutrino scenario one gets 



< 9 eV. In this case we derive from Eq. (|4j) the following branching 

ratio 



Kuu = ^^JT^ ^ f ; < 3.0 x 10~ 30 (3 light neutrino scenario). (6) 
T(K + — ► all) 



Assuming the existence of heavy (i.e. order ~ GeV) mass neutrinos N, we may 
obtain an upper bound on T(K + — > fi + fi + ir~) using the current LEP limit on 
heavy stable neutral leptons M N > 39.5 GeV |T7J, which leads to KM^ 1 )^ < 



n (39.5 GeV) , where n is the number of heavy neutrinos. This limit being 
substituted in Eq. (Q) results in the upper bound 

T^mj, < 2.0 x 10~ 19 (3 light + 1 heavy neutrino scenario). (7) 



Direct searches for K + — ► ji + n + n decay by E865 experiment at BNL JT3] give 

n m < 3.0 x 10" 9 ( 90%CL, Ref. [15]) (8) 

Comparison of this experimental bound with the theoretical predictions in Eqs. 
(PD) (0) cl ear ly shows that both cases are far from being ever detected. On 
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the other hand experimental observation of K + — > {i + fi + 7r~ decay at larger rates 
would indicate some new physics beyond the SM, or, as we will see, the presence of 
an extra neutrino state Vj with the mass in the hundred MeV domain. Let us note 
that such neutrinos are not excluded by the LEP neutrino counting experiments 
measuring the Z-boson invisible width. These experiments set limits not on the 
number of light massive neutrinos but on the number of active neutrino flavors 
N v = 3. Thus extra massive neutrino states Vj can appear as a result of mixing of 
the three active neutrinos with certain number of sterile neutrinos. These massive 
neutrinos are searched for in many experiments [18]. The Uj states would manifest 
themselves as peaks in differential rates of various processes, and can give rise 
to significant enhancement of the total rate if their masses lie in an appropriate 
region. 



3 Hundred- MeV neutrinos in K + —> /i + /i + 7r -decay. 
Resonant case. 

Assume there exists a massive Majorana neutrino Uj with the mass rrij in the 
range 

y^T « 245 MeV < m, < ^4 w 389 MeV. (9) 

Then the s-channel diagram in Fig. 1(a) blows up because the integrand of the 
first term in Eq. ([!]) has a non-integrable singularity at s — m|. Therefore, 
in this resonant domain the total z/j-neutrino decay width T v j has to be taken 
into account. This can be done by the substitution rrij — > rrij — (i/2)T v j. The 
total decay width T v j of the Majorana neutrino Vj with the mass in the resonant 
domain (§) receives contributions from the following decay modes: 

fe + 7T~, e~7T + , /U + 7T~, /U~7T + , 
e + eru c e , e+fi-js*, /i+e"z/ e c , ^[Tv^ , (10) 
e"e + z/ e , e"yu + ^, /i"e+z/ e , yT^v^. 

Since Vj = z/| it can decay in both Vj — > l~X(AL = 0) and Uj — > l + X c (AL = 2) 
channels. Calculating partial decay rates we obtain 

rfa - In) = m'^firntFfayJ = \U V \ 2 T$>, (11) 
V[y . ^ kkv) = p h .^^3_ m p {yiuyi2) = \U hj \ 2 Tf\ (12) 
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where yi = rrii/mj and 

F(x, y) = A 1/2 (l, x\ y 2 )[{\ + x 2 )(l + x 2 - y 2 ) - Ax\ (13) 

H(x, y) = 12f—(z- y 2 ){l + x 2 - z)X 1/2 (l, z, a^A^O, y 2 , z). (14) 



z 



The integration limits are Z\ = yf , z 2 = (1 — y/J 2 and F(0, 0) = H(0,0) = 1. 
Summing up all the decay modes in (^) one gets for the total z/j width 

v vj = 2|[/ w f(r^ + r^ e) + rf M) ) + 2|/y e /(rf } + rf e) + r^) = 

= \U M \ 2 Y^ + \U e3 \ 2 Y^. (15) 



In the resonant domain @ T^- reaches its maximum value at = \j s~{ . As- 
suming for the moment \U^j\ = \U e j\ = 1, we estimate this maximum value to 
be T U j pa 4.7 x 10 -10 MeV. Since T U j is so small in the resonant domain (§) the 
neutrino propagator in the first term of Eq. ([!]) has a very sharp maximum at 
s = m 2 . The second term, being finite in the limit T u j = 0, can be neglected in 
the considered case. Thus, with a good precision we obtain from Eq. (JTJ) 

r -(K+ -> /X+/X+7T-) « C7rG(z ) H^Eml (i 6 ) 

with z = (raj/mx) 2 - This equation allows one to derive, from the experimental 
bound of Eq. (§), the constraints on z/j neutrino mass rrij and the mixing matrix 
elements U^j, U e j in a form of a 3-dimensional exclusion plot. However one may 
reasonably assume that |Z7 W -| ~ l^eil- Then from the experimental bound @ we 
derive a 2-dimensional — |?7 W -| 2 exclusion plot given in Fig. 2. For comparison 
we also present in Fig. 2 the existing bounds taken from ||17|| . As shown in 
the figure, the experimental data on the K + — > /i + yU + 7r~ decay exclude a region 
unrestricted by the other experiments. The constraints can be summarized as 

\U^\ 2 < (5.6 ± 1) x 1(T 9 for 245 Mev < m s < 385 MeV, (17) 

The best limit \U^\ 2 < 4.6 x 1(T 9 is achieved at rrij pa 300 MeV. Note that these 
limits are compatible with our assumption that \U^j\ ~ \U e j\ since in this mass 
domain, typically \U e A 2 < 10~ 9 



The constraints from K + — > /i + yU + 7r~ in Fig. 2 and Eq. (K7h can be signifi- 
cantly improved in the near future by the experiments E949 at BNL and E950 
at FNAL [Tj|. It is important to notice that in the resonant domain we have 
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T res( K + _^ n+^+n-) ~ \u\ 2 , while outside T(K+ -> n+n+TT-) ~ |£/| 4 . Thus in 
the resonant mass domain the K + — > fi + fi + 7i~ decay has a significantly better 
sensitivity to the neutrino mixing matrix element. In forthcoming experiments 
the upper bound on the ratio in Eq. (j8|) can be improved by two orders of mag- 
nitude or even more. Then this experimental bound could be translated to the 
limit |f^ w | 2 < 10~ n and stronger. 

Finally, it is important to mention the possible cosmological and astrophysi- 
cal consequences of the hundred-MeV neutrinos considered in the present paper. 
Massive neutrinos contribute to the mass density of the universe, participate in 
cosmic structure formation, big-bang nucleosynthesis, supernova explosions, im- 
print themselves in the cosmic microwave background etc. (see |20| for a review). 
This implies certain constrains on the neutrino masses and mixings. Currently, 
for massive neutrinos in the mass region of Eq. @ , the only available cosmological 
constraints arise from the mass density of the universe, r v . < (~ 10 14 ) sec, and 
cosmic structure formation, T Uj < (~ 10 7 ) sec, (taken from ref. f20|| ) where r v . is 
the lifetime of massive neutrinos. On the other hand, on the basis of the formula 
flT5|), assuming \U^j\ 2 ~ \U e j\ 2 < 4.6 x 10" 9 as in Eq. fllTf), we find conservatively 
10~ 2 sec < t Vj . Thus massive neutrinos with masses in the interval (|[J) do not con- 
tradict the known cosmological constraints and there remains a wide open interval 
of allowed mixing matrix elements: (~ 10 -18 ) < \U ej \ 2 < (~ 10 -9 ). Big- 

bang nucleosynthesis and the SN 1987A neutrino signal may lead to much more 
restrictive constraints. Unfortunately, as yet the analysis |21| of these constraints 
does not involve the mass region (|9|). It may happen that these constraints, in 
combination with our constrains in Eq. (|ITD , close the window for neutrinos with 
masses in the interval (H). Then the only physics left to be studied using the 
K + — > /i + yU + 7r~ searches would be physics beyond the SM other than neutrino 
issues. Nevertheless, significant model dependence of all the cosmological con- 
straints should be carefully considered before such a determining conclusion is 
finally drawn. 



4 Conclusion 

We studied the potential of the K-meson neutrinoless double muon decay as a 
probe of the Majorana neutrino masses and mixings. We found that this process 
is very sensitive to the handred-MeV neutrinos Vj in the resonant mass range (^j. 
We analyzed the contribution of these neutrinos to the K + — > fi + fi + 7r~ decay 
rate and derived stringent upper limits on the Majorana neutrino mixing ma- 
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trix element |i7 w -| 2 from the current experimental data. In Fig. 2 we presented 
these limits in the form of a 2-dimensional exclusion plot, and compared them 
with existing limits. The K + — > fi + fi + ir~ decay excludes a domain previously 
unrestricted experimentally. We point out that the current and near future ex- 
perimental searches for this decay are not able to provide any information on the 
light eV-mass (m„ ~ eV) or heavy GeV-mass (m v ~ GeV) neutrinos, since in 
those cases the required experimental sensitivities are by many orders of magni- 
tude far from the realistic ones. Finally, we notice that the decay K + — > {i + fi + 7i~ 
can in principle probe lepton number violating interactions beyond the standard 
model. A well known example is given by R-parity violating interactions in su- 
persymmetric models. However these aspects of the K + — > fi + fi + 7T~ decay are 
yet to be studied. 
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(a) (b) 

Figure 1: The lowest order diagrams contributing to K + — ► {i + fi + 7i~ decay. 
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Figure 2: Exclusion plots in the plane \U^j\ 2 —irtj. Here U^j and rrij are the heavy 
neutrino Uj mixing matrix element to and its mass respectively. Domains above 
the curves are excluded by various experiments according to the recent update 



in Ref. |[L7|| . Region excluded by K + — > fi + ^ + 7c decay [present result] covers the 
interval 249MeV< rrij < 385MeV and extends down to \U^\ 2 < 4.6 x KT 9 . 
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